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Abstract. In this article, an uniform discretization of stochastic in- 
tegrals /p 1 f'_(Bt)dB t , with respect to fractional Brownian motion with 
Hurst parameter H £ (5, 1), for a large class of convex functions / is 
considered. In [[I], Statistics & Decisions, 27, 129-143] , for any convex 
function /, the almost sure convergence of uniform discretization to such 
stochastic integral is proved. Here we prove L r - convergence of uniform 
discretization to stochastic integral. In addition, we obtain a rate of 
convergence. It turns out that the rate of convergence can be brought 
as closely as possible to H — 4. 
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1. Introduction 

It is well-known that fractional Brownian motion B = {Bt}te[o,i] ; with Hurst 
parameter H 7^ ^ is neither a semimartingale nor a Markov process. There- 
fore, according to Bichteler-Dellacherie theorem, the classical Ito stochastic 
integration theory cannot be used to define a stochastic integral with respect 
to fractional Brownian motion. In last decade, many authors studied dif- 
ferent possible ways to define stochastic integrals with respect to fractional 
Brownian motion. Essentially two different types of integrals can be defined: 

• The pathwise Riemann-Stieltjes integral L utdBt exists if the inte- 
grand stochastic process u = {ut\t&\o,i\ h as Holder continuous sam- 
ple paths of order a > 1 — H, as a result of Young integration theory 

. Zahle [12] extended this integral using the fractional integration 
by part formula, to some stochastic processes having some fractional 
smoothness. 

• The Skorokhod integral ( or divergence integral )with respect to frac- 
tional Brownian motion. This integral is defined as adjoint operator 
of the Malliavin derivative. It is known that for enough regular sto- 
chastic processes, the difference of this integral with corresponding 
pathwise Riemann-Stieltjes integral can be explained with Malliavin 
trace operator. For more details see [7]. 

In [2] , the authors studied the problem of finding a discrete approximation 
of the stochastic integral with respect to fractional Brownian motion defined 
as divergence. They provide a discrete approximation of stochastic integrals 
of divergence type by means of the resolutions of the Fock space associated 
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to fractional Brownian motion B. 

The aim of this paper is to study in more details convergence of the uniform 
discretization of pathwise stochastic integrals 

(1-1) S= f fUBt)dB t , 



Jo 

where / : 1R — > IR is a convex function. In [TJ, it is shown that such 
integrals can be understood in the generalized Lebesgue-Stieltjes integral 
sense. Moreover, the authors considered the uniform discretization 

n 

S n :=y2fUBi = i)(B i -B i= i), n G IN. 

' " n n n 

i=l 

They proved that S n converges to stochastic integral 5 almost surely as 
n tends to infinity. Note that, from financial application point of view, it 
is convenient to take the left hand points in the definition of S n . In 
stochastic finance the discretization S n can be interpreted as total losses or 
gains of the discretized delta hedging strategy (see [1] ) . 

In this paper, we show that with some fine and detailed analysis one can 
prove U- convergence too, i.e. 

S n — > 5 in L r , 

as n tends to infinity for some range of r > 1 and a large class of con- 
vex functions. To obtain such result, we use integration theory known as 
the generalized Lebesgue-Stieltjes integration theory, introduced by Zahle 
[12] , and developed by Nualart-Rascanu in [8] together with Lemma 13.21 
Moreover, we obtain a rate of convergence: 



, H- 

\r 1 



:iL|.V„ - 51' ] ' <C Q 1 n > 1, 

for sufficiently small e and C is a constant independent of n. 

The paper is organized as follows. In the section 2, we state our main result. 
Section 3 contains all auxiliary facts which we need to prove our main result. 
The section 4 is devoted to the proofs. 

2. Main result 

Throughout the paper, B = {-Bt} tG [ ,i] stands for a fractional Brownian 
motion on the interval [0,1] with Hurst parameter H £ (f>l)- Let / : 
IR — > IR be a convex function, and denote by fi the positive Radon measure 
corresponding to its second derivative. Let p be a positive number such that 

(2.1) 2H<p <Y ^. 
Define a function C : IR — > IR by 

min{a 2 ,(a-l) 2 > 

(2.2) C{a) = max(l, \a\)e ? . 
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We consider the following technical assumptions related to the measure \x: 

(Hi) For p and C(o) given in O and it holds 

i 

C(a)p fj,(da) < oo. 



(H2) For p and C(a) given in 12.11 and 12.21 it holds 



oc 



C(loga)p/u(da) < 00. 







In what follows, the stochastic integrals with respect to fractional Brownian 
motion are understood in the sense of the generalized Lebesgue-Stieltjes in- 
tegrals (see section 3 for more details). Now we can state our main theorem. 

Theorem 2.1. Let f : JR — > 1R be a convex function. Put 

n „i 

Sn = J2fUBi=i){B ± -Bi=i), and S = f_(B u )dB u . 

* ■* n n n f n 

i=l Ji} 
Let r G [1,2?), for p given in h2.1\) . If the assumption (Hi) holds, then 

(a) : We have the convergence S n — > S in L r . 

(b) : For every positive number (3 satisfying 

(2.3) l-F</3<— , /3/1 , 



there exists a constant C = C(H, f3,r,p, f) such that 
\\S n -S\\ r = (JB\S n -S\ r y < C I -J , n>l. 



Remark 2.1. For r G [l,p), the assumption (Hi) implies that S,S n G F r . 

Corollary 2.1. Let X = {^t}tg[o,i] ^ e a geometric fractional Brownian 
motion, i.e. Xt = e Bt . Put 

n „! 

S n = Y J f-{Xi= L )(,X ± -Xi= L ), and S= f'_{X u )dX u . 
i=l ,/u 

Fei p and /3 6e positive numbers such that \2.1\) and \2. 3\) holds and let 
r G If the assumption (H2) holds, then 

(a) : VFe /ioue i/ie convergence S n — > 5 in F r . 

(b) : There exists a constant C = C(H, j3,r,p, f) such that 

||5„-5|| r = (E|5 n -5r) 7 <C(-) > ™> L 
Remark 2.2. For r G [l,p), i/ie assumption (H2) implies that S,S n G F r . 
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Remark 2.3. According to Theorem 2.1, rate of convergence can be any 
number in a certain interval. Therefore, one can improve the rate by choos- 
ing suitable parameters p and /3. In the best, rate of convergence can be 
brought as closely as possible to H — |, by letting p very close to 2H and 
(3 very close to 1 — H . However, as a price the constants C and C become 
larger and tends to infinity as p decreases to 2H or f3 decreases to 1 — H . 

Remark 2.4. It is not clear whether the approximate rate H — \ is the best 
possible rate that one can obtain. 

2.1. Related results and comparison with Brownian motion. For 

more smooth integrands one can use Young-Loeve estimate to get a bet- 
ter rate. In following theorem, the stochastic integral coincides with the 
Riemann-Stieltjes integral (see remark [ 



Theorem 2.2. Let f : IR — > 1R be a Lipschitz function. Put 

n -i 

S n = Y / f{B i=1 ){B±-Bi =1 ), and S= f(B u )dB u . 

Let r > 1. Then for every e € (0,2H — 1), there exists a constant C 
C(e, H,r, f) such that 

\\S n -S\\ r = (TE\S n -S\ r y <C - , n>l. 



For standard Brownian motion W = {Wt} t £[o,i} i«e. H = ^, the stochastic 
integral is understood in the sense of Ito integral. In this case, we consider 
the following assumption related to the measure \x: 



(H 3 ) It holds 

min{a 2 ,(q-l) 2 } 

e 2 /i(da) < oo. 

'R 

Then, we obtain the following result: 



Theorem 2.3. Let W = {Wt}te[o,i] be a standard Brownian motion and 
f : IR — t- IR be a convex function. Put 

n „! 

T n = J2f-(W i= i)(W ± -Wi =1 ), and T= f'_(W u )dW u . 
Let r S [1,2]. If the assumption (H3) holds, then 



(a) : We have the convergence T n — )■ T in L r . 

(b) : There exists a constant C = C(r, /) such that 

\\T n - T\\ r = (JE\T n - T\ r ) r<c(^\\ n > 1. 

Remark 2.5. The rate n~i obtained in Theorem \2.3\ is sharp for the class 
of convex functions satisfying the assumption (H3). For example for the 
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convex function f(x) = (x — a) + , there exists a constant C = C(f) such that 
we have 

C n~* < (]E\T n -T| 2 )^, n > 2, 
where T n and T are as in Theorem WIl 



Remark 2.6. To compare with the case of fractional Brownian motion, we 
see that the rate of convergence is better for fractional Brownian motion 
than for standard Brownian motion if H > |, and worse if H < |. It is 
known that (see [3]) the mixed Brownian-fractional Brownian motion X = 
W + B is a semimartingale, if H E (|, 1), and for H £ | ] , X is not a 
semimartingale with respect to its own filtration IF . 

For smooth functions, using ltd isometry one can easily get the rate n~2. 
This is the subject of the next theorem. It is known that this is the best 
possible rate one can get for Lipschitz functions. See [TQl Remark 3, p. 694] 
and references therein. 

Theorem 2.4. Let W = {Wt} t £[o,i\ be a standard Brownian motion and 
f : JR — )■ IR be a Lipschitz function. Put 

n „i 
T n = Y / f(Wi =1 )(W ± -Wi =1 ), and f= f(W u )dW u . 
i=i J{J 

Let r > 1. Then there exists a constant C = C(f,r) such that 

\\f n -f\\ r = (m\f n -f\ r )r < C^=, n>l. 

3. Auxiliary facts 

3.1. Pathwise stochastic integration in fractional Besov-type spaces. 

Since fractional Brownian motion in not a semimartingale, hence the sto- 
chastic integral with respect to fractional Brownian motion B must be de- 
fined. Using the smoothness of the sample paths of B, when H £ (^, 1), one 
can define the so-called generalized Lebesgue-Stieltjes integral. We shall give 
some details of the construction of generalized Lebesgue-Stieltjes integrals 
in this section. For more information see [3 Section 2.1.2]. 

Definition 3.1. Fix < /3 < 1. 

(i) Let = W^QOjT]) be the space of real-valued measurable functions 
f : [0, T] ->• R such that 

\\f \\ip := sup — 7g—+ — ( n+B~ du <°°- 

(ii) Let W2 = W2 ([0, T]) be the space of real-valued measurable functions 
f : [0, T] -)■ M such that 
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Remark 3.1. The Besov spaces are closely related to the spaces of Holder 
continuous functions. More precisely, for any 0<e</3A(l — /3), 

C^ +e ([0,T]) C Wf ([0,T]) C C^([0,T]) and C^+ e ([0,T]) C Wf([0,T]) ; 
where C 7 ([0,T]) denotes Holder continuous functions of order 7. 

Recall that the trajectories of B belong to C 7 ([0,T]) almost surely for any 
T > and any < 7 < H. This follows from the Kolmogorov continuity 
theorem. By remark 13.11 we obtain that the trajectories of B belong to 
Wi ([0,T]) almost surely for any T > and any < j3 < H. 

Denote by T the Gamma-function. Recall the left-sided Riemann-Liouville 
fractional integral operator of order (5 > 0: 

(i&/)W = ^ ) [f(u)(s-uf- l du. 
The corresponding right-sided fractional integral operator I_ is defined by 

(llf)(s) = ^ f(u)(u- sf- l du. 

Remark 3.2. Iff G Wf ([0, T]), i/iera ife restriction to [0,t] C [0,T] 6eZon#s 
to /^(LooftO,*])). Afeo, i// G W|([0,T]), £/ien its restriction to [0,t] C [0, T] 

6e/on 5 s to /^(^([O,*])), w/iere /^(^([O, t])) (Vesp. (Li([0, stand for 
the image o/L oo ([0,t]) (resp. L\{[0,t})) by the fractional Riemann-Liouville 
operator I_ (resp. I\L).(For details we refer to [SJJ. 

Definition 3.2. Let f : [0, T] ->■ R and < /3 < 1. /// G jJ(Li([0, T])) (resp. 

/ G /^(Loo([0, T])) i/ten i/ie Riemann-Liouville fractional derivatives are de- 
fined using the Weyl representation as 



(rfi *\( \ 1 ffW 1 r f x f( x )-f(y) , \ 1 f, 

iD °+ f){x) = W^W)\^ r + PJo {x-y)W dV ) 

[resp.{D T _f){ X ) = ^-^ [ W ^ W + Pj x iy _ x)P+ i d V) W*)J " 

For a detailed discussion, we refer to [Jj]. The following proposition clarifies 
the construction of the stochastic integrals. This approach is by Nualart 
and Rascanu. 

Proposition 3.1. [8j Let f G W$([0,T\), g G W^([0,T]). Then for any 
t G (0, T] the Lehesgue integral 

!l(D^ + f){x)(D l - p g t -)(x)dx 
exists, and we can define the generalized Lebesgue-Stieltjes integral by 



fdg:= / (D^ + f)(x)(D^ p g^)(x)dx. 



RATE OF CONVERGENCE 7 

Remark 3.3. It is shown in [12j that if f G C^([0,T]) and g G C^([0,T]) 
with'y + fi > 1, then the integral fdg exists in the sense of the generalized 
Lebesgue-Stieltjes integral and coincides with the Riemann-Stieltjes integral. 

The next theorem can be used to study the continuity of the integral. 

Theorem 3.1. [5] Let f G W$[0,T\ and g G W^^[0,T}. Then we have the 
estimation 

rt 



(3.1) 



fdg 



< sup I^PVOO \\\fh,p, te[0,T). 

0<s<t<l 



Now we can state the existence of stochastic integral with respect to frac- 
tional Brownian motion in our main result. 

Theorem 3.2. pQ Let f : JR — )■ TR be any convex function, 
(i) The stochastic integral 

f'-(B t )dB t 



o 



can be understood a.s. in the sense of the generalized Lebesgue-Stieltjes in- 
tegral. 

(ii) The following Ito formula 

/(B0 = /(0)+ / f'-(Bt)dB t 
Jo 

holds, where the stochastic integral is understood in the sense of the gener- 
alized Lebesgue-Stieltjes integral. 

(Hi) One can approximate the stochastic integral by Riemann-Stieltjes sums. 
More precisely, 



n .% 

J2f'-(B t ? JiBt* - B t? ) ^ / fUBt)dB t , tf 

Jo 



i 

n 



i=l 

3.2. Some results related to fractional Brownian motion. The so- 
called Garsia-Rademich-Rumsey inequality provides basic inequalities on 
increments of continuous stochastic processes. Using this inequality, one 
can obtain the following lemma on the moments of supremum of fractional 
derivative of fractional Brownian motion. 

Lemma 3.1. [8] Let B = {Bt}t^\o,i] be a fractional Brownian motion with 
Hurst parameter H G (5, 1). Let 1 — H < (3 < \ and p > 1, then 

E( sup \D}r P B t -(s)\ Y < 00. 

0<s<t<l 

We continue with an useful estimate of a probability that fractional Brow- 
nian motion crosses a fixed level. It turns out that this is a main ingredient 
for the proof of the main theorem. Actually the following result is an im- 
provement of the Lemma 4 (see [1]) with a better constant in terms of the 
level a. 
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Lemma 3.2. Let B = {Bt}te[o,i} be a fractional Brownian motion with 
Hurst parameter H G (|,l)- Fix < s < t < 1 and a G H. T/ien 
i/tere exists a constant C , independent of s, t and a, such that the following 
estimate 

TP(B t > a and B s < a) < C C{a){t - s) H s~ 2H 

holds. 

Lemma 3.3. Let W = \Wt\te\a,i] be a standard Brownian motion. Fix 
< s < t < 1 and a G IR. Then there exists a constant C, independent of 
s, t and a, such that the following estimate 



, s mm{a 2 ,(g-l) 2 } If 

JP(W t > a and W s < a) < Ce * 

holds. 

The proof of the lemmas are given in Appendix A. We also use the follow- 
ing well-known estimate for the tail probability of standard normal random 
variable. 

Lemma 3.4. Let Z be a standard normal random variable and fix a > 0. 
Then 

(3.2) W{Z>a) < -^e"f 



4. Proofs 

We start with the following simple lemma. It turns out that it provides 
enough good upper bound. 



Lemma 4.1. Let n > 2 and a G (0, 1). Then 



i=l 



—n 

a 



Proof of theorem \2.1\ Throughout the proof all constants will be denoted by 
C, and their values may differ from line to line. Random constants will be 
denoted by C(w). We prove the statement only for r = 1. The general case 
follows by similar arguments (see Remark l4.1|) . Note that 

Sn~S= [ h n (t)dB t 

Jo 

where 

n 

(4.1) h n {t) = £ (/i(B^i) - fUBtj) l(i=l,i](<)- 

i=l 

By Theorem 13.11 and Lemma 13.11 there exists a random variable C{uj,H,(3) 
for which all the moments exists and 

(4.2) \S n -S\<C{u,H,0)\\h n h,p 

for every /9g(1 — H, h). Thus by Holder inequality, we obtain 



(4.3) JE\S n -S\< C(H,p,p) JE\\h n \\ p 2 
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Let now p be as in (|2,ip and let j3 £ ( 1 — H, ^ ) . We proceed to compute 



the term 



mhn\\ P 2d 



P . We have 



<QEJV)p + QElP)t, 



where J n denotes the first term and I n the second term in the Besov norm 
|| • \\28- The rest of the proof is split into three steps. We first prove the 
statement for a convex function f(x) = (x — a) + , where o £ 1R. Next 
we prove the statement for convex functions for which the measure \i has 
compact support. Finally, we prove the result for convex functions for which 
the assumption (Hi) holds. 

Step 1. The case fix) = [x — a) + . 

Now we have 



i=l 



(4.4) K(t) = ( l{ J B t <a<B i _i} - l{B i _ i <a<B i } ) ,i] 



(*)• 



For the term J n , we use Minkowski inequality for integrals to obtain 
- 1 (JE\h%(t)\P)p 



M)5 < 



< 



W{B t > \a\) l 'P 



dt 



» ri TP(B t >a> B^y/P + W(B t <a< B^i) 1 ^ 

+ y: / * — -» = — ^ 



i=2 

:= JnA + Jn,2- 



IP 



In J n> i, the probability can be estimated by one, if |a| < 1, and by estimate 
3.2 if |a| > 1. Hence the term J n i can be bounded as 



1 



Jn l < C / -jdt < C - 



o f 



1-/3 



For the term Jn^j by symmetric property of fractional Brownian motion, it 
is sufficient to consider only the event {B t < a < B^}. Therefore, J n ^ can 

n 

be bounded using Lemma 13.21 as 

Jn,2 < / i ^ n) .} n) dt 



i=2 



^ W V n 
/ i \ 1-—-P n 



-2H/p 



n 1 

-~dt 

»-i t/ 3 
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Together with lemma B~Tl this implies that 



j si C &)* if " €(0,1), 

n ' 2 - 

C(i) ifa>l, 



where a = ^ + /3. 

We proceed to study the term I n . We split the integral into several parts. 
Particularly, we consider the cases when s and t lie in the same interval 



«, Ci VCUO. .,wu C uxxcu VVXX^ o c y 

* G £ I with * ^ J> we have 



and when they lie in different intervals. Note that when s E ( ^— , ^ 



and 



K(t) -/£(«)! 

= l 1 {Bi<a<Bi_i } - l{B t >a>B i _i} ~ 1{B S <a<Bj_i } + l{B s >a>Bj_i } I 

n n n n 

< l 1 {B i <a<B 1 _i } - l{B s <a<B 2 _i }l + 1 1{B S >a>B i _i } ~ l{B t >a>.Bi-i } I 

n n n n 

:=H 1 (j,i)+H 2 (j,i). 
Using Minkowski inequality for integrals, we have 



JO 



. r* P(B t > a > B s y/P + P(B t < a < B s ) 1 l p A+ 

- ' 1 (T^W 1 



J 



» fi [> P(B t > B >B.)-/P + P(B,<,<B.)V 
^ /.-i (t - s)^ 1 

Q)j + (11 

(t - s)^ 



i=2 



n n 



(t - a)* 



n ** n 



i=2 j=l 
In,l + In,2 + ^n,3 



We start with I n ^. Note that it is enough to consider only the term Hi(j, i). 
The term H 2 (j, i) can be treated similarly. We have 



i / \ lip / \ l/p 

(JEHf(j,i))p < 2P (B s < a < Bj-i) + 2P ( B t < a < Bi-i J 
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Hence the term 3 can be bounded as 



n i—1 ,.L ,.L 
i=2 j=l 



IP ! D„ < « < D, , ) 1/P + P ( B t < n <- B 



i/p 



_ 3=2. 

7i " n 



n i—1 



c ee 



i=2 j=l 



n n 



(<-*) 



3+1 



(t - s)/^ 1 
1/p 

— dsdt 



-dsdt 



\i/p 



+ Ji=i J i=i ft - 5)^+1" 

i=2 j = l ^ — n V y 

- + / (2) 

For the term J v 3, by using Lemma 13.21 we obtain 



-dsdt 



s^e/,: 



i=2 
n 



i-1 
n In 



2H 
p 



/9+1 



< 



c e;, 



i=2 



(t-s) 



i — 1\ p 



-dsdt 



n 



n 



2H 
V 



i - 1 



11 



dt 



n 



where for the last inequality, we have used Lemma 14.11 Next we consider 
the term In this case, we have to study the case j = 1 separately. Let 

j = 1 in the term i^g. Then, by proceeding as for J n> ±, we have 



Ef, 



i=2 



1/p 



» P (g, < a < 0) 
(t - s)^ 1 

n 

^e /: 



■dsdt 



< 



i=2 
-1 



rt I n 



(t-s) 



J+i 



dsdt 



1 

t 

n 



dt 



71 



1+1 - 

n 



-C|i 

n 



1- 



If j > 1, then by changing the order of two summations and tedious manip- 
ulation, one gets 
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— dsdt < C 



eel; 



it - sf 



+1 



. n 

i=2 j=2 ~ " 

It remains to estimate the terms I n) \ and l n ^i- For I n i, by using Lemma 
we have 

(t — s)p s p 



= cy tT J s (}~t) dsdt 

= C/ (tu)~T (1 -n)^^ _1 dndt 

jo •/ o 

2H H fn H a 

= C B{1 - — , — - p) / r p _/, dt 
P P Jo 



C 



n 



where B(x,y) denotes the complete Beta function. For the term I n 2, we 
obtain 



' 1 (t — S) p S P , , 

dsdt 



" r— ft 



n n 



i=2 
n 



' Z — 1 \ P n I I — 1\ P 



i=2 

cfi 



/ 1 \ 1— — — j8 n / i \ — 



<C|- 

71 



1 \ P 



where we have used Lemma B~T1 Finally, by collecting estimates for J n ,i, 
Jn,2, In,!-, I-n.,2 and I n> 3 we obtain that for the convex function f(x) = (x — 
a) + , there exists a constant C = C(H,f3,p) such that 

1 (\\-* < p 

(4.5) JE\S n -S\<C C(a)p [- 

\n 

Step 2. The case supp(//) is compact. 

It is well-known that left derivative of the convex function / has the following 
representation 

f'-( x ) = n I s S n (. x - a)fi(da) 
1 Jn 
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up to constant. So 
1 n 

i=l 
1 n f 

= 2 E / ( s S n (- B ^i ~ a ) - sgn(5 t - a)) //(da)l^ ^ (t) 
i=r E 

1 f n 

2 je ~[ L -r -r 

= [ h a n (t)n(da). 
Jr 

From this observation, we obtain 
(4.6) 
and 

(4.7) \h n (t) - h n (s)\ < I \h a n (t) - h a JsMda) 



\ n ' n J 



\K(t)\ < / |^(t)| M (da) 
JE 



Hence, using Tonelli's theorem, Minkowski inequality for integrals and in- 
equality (|4.6p we obtain 



and 



IE 



IE 



< 



R JO 



1 (E|/j « (t)|P) "-dtMda), 



1 r t 



dsdt 



o jo 

< 



\K(t) - K(s)\ 

(t - s)P +1 

1 *QE\hZ(t)-htt8)\P)i 



dsd*/i(do). 

RJO JO I 1 - s r 

Therefore, using step 1 we can conclude that there exists a constant C = 
C(H,/3,p) such that 

— — 8 

E\S n -S\ <C J C(a)*iM(da) f~V , 

where C(o) is given by (|'2.'2p . 
Step 3. The general case. 
Now take any convex function / which satisfies the assumption (Hi). For 
any k £ IN, define the measurable set by 



(4.8) 



fifc = {uj : sup 6 [0, k]} 

0<t<l 



and auxiliary convex functions by 

'fU-k)x + (f(-k) + fL(-k)k), x<-k 
f(x), x e 

U+(*)* + (/(*0-/+(*0*0> 



(4.9) / fc (x) = ^ 
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Denote by the positive Radon measure associated to the second derivative 
of convex function Then /z^ has compact support contained in [—k,k]. 
Let S k and S k stand for the stochastic integral and the uniform discretization 
as in main theorem corresponds to the convex function fy. Note that on the 
set fifc, we have S k = S n and S k = S almost surely. Hence, by monotone 
convergence theorem, we have 

E|5„-5| = lim JE\S n -S\l Qk 

= lim JE\S*-S k \ln k 

< lim JE\S k - S k \. 

k— >oo 

Applying step 2 and assumption (Hi) completes the proof. □ 

Remark 4.1. The result for 1 < r < p follows with the same argument by 
choosing suitable parameters in Holder inequality. 

Proof of corollary \2.1\ We prove the result for the function f(x) = (x — a) + 
with some positive constant a. For negative a, we have (Xt — a) + = Xt — a 
and the result is trivial. Moreover, the result for general convex function / 
satisfying the assumption (H2) follows by same arguments as in the proof 
of Theorem Put 

(4.10) h*> a (t) = ( 1 {^>a} " **l(i=l f i](t). 

1=1 \ " / 

Then it follows that 

S n -S= f h*' a (t)dB t . 

Jo 

A simple calculation gives us 

hX> a (t) = X t h l °* a (t) 
where h^{t) is given bv 14.41 Hence for the first term J n , we obtain 

■i^wi dt£X /"rt d , 



«« - h * 

where X = sup <i<i Xp Moreover, all moments of X are finite (see [6]). So 
we can replace C(uj,H, f3) in the inequality (|4.2p by a new random variable 
C(u,H,(3) = C(u, H, j3)X . Hence the result follows by Step 1 of the proof 
of Theorem 12.11 Next we consider the second term I n . Note that 

\h*> a (t) - h*> a (s)\ < \ti°* a (t)\\X t - X s \ + \X s \\h^ a (t) - C ga (*)|- 

For the term \X s \\hn ga (t) - h l ° ga (t)\, we can proceed as for J n in the proof 
of Theorem EU For the term \h l ° ga {t)\ \X t - X s \, for any G (0,H - (3), 



JO 
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using the Holder continuity property of sample paths of Xt, we obtain 

(t - s)W d3dt 

<C(co) [ \ti° sa (t)\ [ {t-s) H -P-P- x 6s6A 
Jo Jo 

where C(oj) = C(oj, H, 0) is a positive random variable for which all mo- 
ments are finite. Hence the result follows from Theorem I2JJ □ 



Proof of theorem \2.SX Since / is a Lipschitz function, there exists an univer- 
sal constant L > such that 

\f(y)- f(x)\<L\y-x\, Vx,yGlR. 

It is also known that sample paths of fractional Brownian motion are of 
bounded p— variation almost surely for any p > hL So, let p, q > jj, and 
using Young-Loeve estimate (see [5]), we have for every < e < H 

n 

\S n -S\<CJ2 ll/(S)|| 8 _ w[4=i>1] ||B|| p _ VBrI i =lf±] 

i=l 

n / -i \ 2H—e 



< cm I i 



i=l 

2H-l-e 



n 



for some positive random variable C(u) = C(uj,e,H, f) for which all the 
moments are finite. Now the claim follows. □ 

Proof of theorem \2.3l The result follows by considering the convex function 
f{x) = (x — a) + and applying ltd isometry and lemma [331 The general case 
follows by the same arguments as in the proof of Theorem 12.11 together with 
the assumption (H3). The details are left to the reader. □ 



Appendix A. Proofs of lemmas 13.21 and 13.31 
We begin with the following simple lemma which we use in the proof. 
Lemma A.l. Let H > \ and fix < s < t < \. Put 

R{t ,s)= 1 -[t 2H + s 2H -(t 
Then there exists a constant C such that 



s 



R(t,s) ~ y ' 
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Proof. Note that since H > ^, we have R(s,s) < R(t,s). Let now i > 2s. 
Then 



s H ~ R(t,s)' 
Hence it is sufficient to consider the case s < t < 2s. In this case we have 

R(s, s) ~ R(s, s) 
Hence we only have to prove that 

i R(s,s) < c (t-s) H R(t,s) 



R(t,s) ~ s~ H R(s,s)' 
By putting k = - and dividing with s 5H , this is equivalent to 

[k 2H — l — (k— l) 2H ] < C{k - if [k 2H + 1 - (k - 1) 2H ] 2 
Now we have k G [1,2]. Hence 

[ k 2H _ 1 _ (fc _ j)2g] 

(fc - l) H [k 2H + 1 - {k - 1) 2H } 2 

k 2H _ 1 k 2H_ 1 
< t^t < 



(k - 1) H ~ k H - 1 
= k H + 1 < 2 H + 1. 

This completes the proof. □ 

Proof of lemma \3.2\ Let s) denotes the covariance function of fractional 
Brownian motion given by 

R{t,s) = \[t 2H + s 2H -{t-s) 2H ]. 
We make use of decomposition 

R(s, s) 

where Y is iV(0, 1) random variable independent of B s and 

2 _ R(t,t)R(s,s) - R(t,s) 2 
° ~ R(s,s) ' 

Assume that 

—7 Ma - 1 < a. 

R(t,s v ; 
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Then we obtain 
JP(B t > a> B s 



P \ Y > 



_ R{t,s) 
a R(s,s) X 



a 



2tts h 



P \Y > 



_ R(t,s) 
a R(s,s) X 

a 



e 2 3 2« dx 



1 , 

e 2 S 2H ^ x 



+ 



p y > 



_ R{t,s) 
a R{s,s) X 



2tts h 
1 



2vrs H 



e 2s 2H dx 



We begin with I\. By lemma l37il we have 



JP\Y> 



_ R{t,s) 
a R(s,s) X 

a 



< 



1 



2irA(x) 



e 2 



where ^4(x) 



it(t,s) 

" -R(3,3) X 



. Hence 



h < 



i 



A(z) 2 1 

e 2 — — e 2 3 2-w d x 



< p 2 

— e H 



x I a 



2-kA[x) V2tts h 

1 

e 2 ^ 2H 2 dx 



2ir 



R(s,s) a _ s 2 f°° 1 _j/ 2 L«( 



e 2 



s) s H J 1 2tt 



-e 2a^ 



2H 



si 

2 dy 



Note that a < (t — s) H and R(s, s) < s), so it remains to show that the 
integral is bounded by a constant independent of s, t and a. It is easy to 
see that 



2a 2 



1 

2^ 



y - a 



R(s,s) 
R(t,sf 



a 2 



-R(S,S) _ 2 _ 2 -R(S,«) 2 _ 4 



where 



Now 



and 



i?(t,s) 



I 1 + R(s,s) 



:a — a 



R(t,s) 



-a 



a 2 a 2 R(t,s) 2 ' 
1 



> 1 



(7 



R(S,S) _ 2 _ 2 i?(s,s) 2 _,; 



R{t,s) 



:(T — a 



R{t,s) 



a* > 0. 
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Hence 



y 2 [ gfefyC-^l , a 2 

il2"H 



< 



< 



1 _y" 

— e 2^ 51273 I" 2 Ay 

00 ^ 



2?r 



1 f R(s,s) _ 2 N 

.g 2^V fl(t,s 2 j 



dy 



2vr 



Hence for ii, there exists a constant C such that 

We proceed to study the term I2. Note that a 2 > 0. Hence 

R(s,s) > 1 >L 



As a consequence, there exists a constant C such that 



e < c e - 



n{q 2 ,(q-l) 2 } 



for every a and every a; G 
/ 2 



fi(s,s) 



(a — 1), a 



Hence 



^(z) V27T 



1 _Z 1 * 2 

e 2 : — — —e 2° 2H dx 



< 



2tts h 



00 1 

>' s j(a-l) 



-R(t, 



v 

e 2 



dydx. 



By applying Tonelli's theorem, the integral can be written as 



e 2 dydx 



+ 



V2T 

[ K(s,s)J er 



1 _Z 

e 2 



a - — r (a - ay) 



i2(t,s) 



dy 



1 _Z 

e 2 



a — 777 — r( a ~~ 1) 
R(t,s) K y 



dy 



/2vr 

= : ^2,1 + ^2,2- 

For 12,2, by lemma [331 we obtain 

J 2)2 < Cmax(l, \a\)a < 2max(l, \a\){t - s) H . 
For I21, by applying lemma |A. 11 we obtain 

(t-s) H 



hi < Cmaxfl, Id 1 rr 



Note that if 



R{s,s) 
R(t,s) 



(a — 1) > a, 



^We have C = 1 except when a g [0, 1]. In this case C = e 2 . 
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then we proceed as for I\ and obtain the result. This completes the proof. 

□ 

Proof of lemma The result follows by the same arguments as in the 
proof of lemma [3721 together with the fact that for standard Brownian motion 
we have R(s, s) = R(t, s) for s <t. □ 
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